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ABSTRACT 

Supermassive black hole binaries (SMBHBs) in galactic nuclei are thought to 
be a common by-product of major galaxy mergers. We use simple disk models 
for the circumbinary gas and for the binary-disk interaction to follow the orbital 
decay of SMBHBs with a range of total masses (M) and mass ratios (g), through 
physically distinct regions of the disk, until gravitational waves (GWs) take over 
their evolution. Prior to the GW-driven phase, the viscous decay is generically 
in the stalled "secondary-dominated" regime. SMBHBs spend a non-negligible 
fraction of a fiducial time of 10^ years at orbital periods between days < torb ^ 
year, and we argue that they may be sufficiently common to be detectable, pro- 
vided they are luminous during these stages. A dedicated optical or X-ray survey 
could identify coalescing SMBHBs statistically, as a population of periodically 
variable quasars, whose abundance obeys the scaling A^var oc t"^^. within a range 
of periods around tvar ~ tens of weeks. SMBHBs with M < 10^ M©, with 
0-5 < « < 1.5, would probe the physics of viscous orbital decay, whereas the 
detection of a population of higher-mass binaries, with a = 8/3, would confirm 
that their decay is driven by GWs. The lowest mass SMBHBs (M < 10^"^ M©) 
enter the GW-driven regime at short orbital periods, when they are already in 
the frequency band of the Laser Interferometric Space Antenna (LISA). While 
viscous processes are negligible in the last few years of coalescence, they could 
reduce the amplitude of any unresolved background due to near-stationary LISA 
sources. We discuss modest constraints on the SMBHB population already avail- 
able from existing data, and the sensitivity and sky coverage requirements for a 
detection in future surveys. SMBHBs may also be identified from velocity shifts 
in their spectra; we discuss the expected abundance of SMBHBs as a function of 
their orbital velocity. 



Subject headings: black hole physics - galaxies: nuclei - gravitational waves 
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Introduction 



Supermassive black holes (SMBHs) app ear to be present in the nucleus of most, and per- 
haps all, nearby galaxies (see, e.g., reviews by lKormendy &: Richstonelll995l and lFerrarese fc Ford 
20051 ). The correlations between the masses of the SMBHs and various global properties of 
the host galaxies suggest that evolution of SMBHs is closely related to the evolution of 
galaxies. In particular, in hierarchical structure formation models, galaxies are built up 
by mergers between l ower-mass progenitors. Each merger event is expected to deli ver the 
nuclear SMBHs (e.g. ISpringel. Di Matteo fc Hernauistll2005l : [Robertson et al.ll2006l ). along 
with a significant amount of gas ( Barnes fc HernquistI 1992 ). to the central regions of the 
new post-merger galaxy. 

There is some evidence for nuclear supermassive black hole binaries (SMBHBs), which 
would b e expected to be pro duced in galaxy mergers. Direct X-ray imaging of an active 



nucleus (iKomossa et al. |2003| ) has revealed a SMBH binary at a separation of ~ Ikpc, and 
Boroson fc Lauerl (120091 ) recently identified a candidate SMBHB, at ~ 10^ times smaller 
separation, from its optical spectrum. A radi o galaxy is also known to have a double 
core with a projected separation of ~ 10 pc (IRodriguez et al.l l2006l ). and several other 



observations of radio galaxies, such as the wiggled shape of jets indicating pr ecession (e.g. 
Roos. Kaastra fc Hummellll993l ). the X-shaped morphologies of radio lobes (e.g. lMerritt fc Ekers 



2002 



[e.g. 



Liu 



200 



the inte r ruption and recurrenc e of activity in double-double radio galaxies 



Schoenmakers 



2000 



Liu. Wu fc Cadl2003l ). and the elliptical motion of the unresolved 



core of 3C66B JSudou et allboOsN I have all been interpreted as indirect evidence for SMBH 



binaries down to sub-pc scales. 

Two interesting conclusions may be inferred from the above observations. First, while 
there is evidence for a handful of nuclear SMBHBs, these objects appear to be rare. This 
suggests that if binaries do form frequently, then they coalesce (or at least their orbital 
separation decays to undetectably small values) in a small fraction of the Hubble time. 
Second, SMBHBs can apparently produce bright emission, with a luminosity comparable to 
active galactic nuclei (AGN), before they coalesce. In general, the circumbinary gas, delivered 
to the nucle us in galactic mergers, can both play a catalyst role in driving rapid SM BHB 
coalescence ( iBegelman. Blandford. &: Reeslll980l : iGould fc Rixll2000l : lEscala et al.ll2004l ). and 
could also accrete onto one or both SMBHs, accounting for bright emission during the orbital 
decay. 



^The lack of any modulation in arrival times for radio pulsars suggests that the elliptical motion of the 
last source has a different origin; see ? 13.31 below. 
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The dense nuclear gas around the BH binary is e xpected to cool rapidly, and set - 
tle into a rotationally supported, circumbinary disk (e.g. 



Barnes 


2002; 


Escala et al. 


2005) 



The dynamical evolution of a SMBH B embedded in such a thin disk has been studied 
in various idealized configurations fe.g. lAr mitage fc Nataraian booj lLiiL_Wu fc Caolbood: 
Milosavlievic fc Phinnevlboosi : IPotti et al.ll2006 : MacFadven fc Milosavlievidboosl : Iflavasaki 



2009 



Cuadra et al.ll2009l ). The generic conclusion of these studies is that initially, the or- 
bital decay is relatively slow, and is dominated by viscous angular momentum exchange with 
the gas disk, whereas at small separations, the decay is much more rapid, and is eventually 
dominated by gravitational wave (GW) emission. 

Whether the decaying SMBHBs produce bright electromagnetic (EM) emission is com- 
paratively much less well understood. If the disk is thin, the torques from the binary cre- 
ate a central cavity, nearly devoid of gas, within a region about twice the orbital separa- 
tion (for a nearly equal-mass binary, e.g. lArtymowicz fc Lubowl Il994j ). or a narrower gap 
arou nd the orbit of the lower-m ass BH in the case of unequal masses q = M2/M1 <^ 1 
(e.g. lArmitage fc Natarajanll2002l ). In the latter case, the lower-mass hole "ushers" the gas 
inward as its orbit decays, producing a prompt and luminous signal during coalescence. In 
the former case, if the central cavity were indeed truly empty, no gas would reach the SMB- 
HBs, and bright emission could n ot be produced. However, numerical simulations suggest 
residual gas infiow into the cavit y ( Artvmowicz fc Lubowl Il996l : iMacFadyen &: Milosavljevic 
20081 : iHayasaki et al.ll2007l . l2008l : ICuadra et al.ll2009l ). which may plausibly accrete onto the 
BHs, producing non-negligible EM emissionj^ Finally, S MBHBs recoil at the time of thei r 
coalescence due to the emission of gravitational waves (iBlanchet. Qusailah fc Willi l2005l ). 
The gas disk will respond promptly (on the local orbital timescale ) to such a kick, which 
may produce shocks, and transien t EM emission, after coal escence ( iLippai. Frei fc Haiman 



2008al : ISchnittman &: KrolikI l2008l : IShields fc Bonnine] 12008 ). Th e kick, however, can begin 



building up during the late inspiral phase (ISchnittman et al. 



20081 ). possibly resulting in some 



emission even before the final coalescence. During the late stages of coales cence, emission 
may also be produced by viscous heating of the disk by the GWs themselves ( iKocsis &: Loeb 
20081 ). 



The luminosity, spectrum, and time-evolution of any EM emission produced by coalesc- 
ing SMBHBs, especially during the last, GW-driven stages, remains uncertain. However, 



^This would be followed by an X-ray "afterglow" ~ 7(1 + z)(M/lO^M0)^-^^ yr after the coalescence, 
caused by the gas outside the cavit y falling in, after a delay set by the disk viscous time ( Liu. Wu fc Gaol 



2OO3I : iMilosavlievic fc Phinnevll2005| ) . Such an afterglow is interesting, for example for a follow-up to SMBH 



merger events detected by the Laser Interferometric Space Antenna (LISA), but not relevant to the idea 
proposed in the present paper. 
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any emission produced during the inspiral stage is likely to be v ariable. For example, recent 



nume rical simulations of an equal-mass binary on parsec scales ( iMacFadven &: Milosavlievic 



2008 ) , and of both equal and unequal-mass binaries on sub-parsec scales ( iHayasaki et al. 



20071 . l2008l : ICuadra et al.ll2009l ) find that the circumbinary gas disk is perturbed into eccen- 
tric orbits by the rotating quadrupole potential of the binary, and that the rate of residual 
accretion across the edge of the cavity is modulated, tracking the orbital period. The lu- 
minosity is likely to be directly tied to the mass accretion rate, and therefore may vary 
periodically. However, even if the gas accretion rate were steady, one would expect periodic 
flux variations, due to the orbital motion of the binary (Kocsis & Loeb, in preparation). 

In this paper, we address the question: Given their expected rate of orbital decay, could 
the population of coalescing SMBHBs be identified statistically in an observational survey for 
periodically variable sources ? Given that the interpretation of individual SMBHB candidates 
have so far remained ambiguous, with alternative explanations possible for each source, the 
potential for such a statistical identification should be explored. 

To answer this question, we first utilize steady-state thin disk models to study the 
orbital decay of a SMBHB, embedded in a circumbinary disk. The decay is described by 
the residence time t^es = —R{dR/dt)~^ the sources spend at each orbital radius R, or at the 
corresponding orbital timescale torb- In the limiting case of a purely GW-driven evolution, 
which becomes valid at small orbital separations (typically at < several x 10^ Schwarzschild 
radii, but with large variations; see below) and remains valid until the final minutes of the 
merger (the so-called "plunge" stage), the residence time is given by tres = ^gw oc t^^. At 
larger separations, the viscous interaction between the binary and the disk drives the binary 
evolution. The residence time in this regime becomes dependent on assumptions about the 
properties of the disk and the nature of the binary-disk interaction, which we will explore in 
this paper. In general, tj-es oc t"^^, with the generic value of a well below 8/3 - significantly 
flatter than the t^^sVS. torb relation in the GW-driven stage. 

We then hypothesize that (i) non-negligible emission (at a fair fraction of the Edding- 
ton luminosity) is maintained throughout the orbital decay, and (ii) the luminosity varies 
periodically on the orbital time-scale. The first assumption allows us to identify coalescing 
SMBHBs with luminous quasars. The second assumption implies that as the orbit of a 
binary decays, its variability timescale decreases. Among sources at redshift z with similar 
inferred BH masses, the observed incidence rate /var of periodic variability on the time-scale 
^var ~ (1 + -2^)^orb, IS then proportioual to the residence time tres = ^res(^var)- At short periods, 
the /var could therefore show a characteristic power-law dependence on tvar indicative of 
a GW-driven evolution, whereas at longer periods (and, as we will discuss, for lower BH 
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masses) the dependence will be flatter, due to viscosity-driven evolutionj^ 

We quantify the requirements that such periodically variable sources be identifiable, 
based on their incidence rate, in an optical or X-ray survey. Luminosity variations at a 
fraction f^dd ^ 0.01 of the Eddington luminosity would correspond to a periodically varying 
flux component with amplitude (/Edd/0.01)(Mbh/3 x 1O'^M0)1O"^^ erg s"^ cm"^ for BHBs 
at 2 = 2, or to i ^ 26 + 2.5 log[(/Edd/0.01)(MBH/3 x IO^Mq)-^] magnitudes in the optical. 
We find that these periodic sources are either too faint or too rare to have been found in 
existing variability surveys. However, if the overall luminosity is indeed a non-negligible 
fraction of the binary's Eddington luminosity, then a long-duration future survey, sensitive 
to periods of weeks to tens of weeks, could look for periodically variable sources, and identify 
a population of sources obeying well-defined scaling laws. 

The discovery of a population of such periodically variable sources could have several 
implications. At long periods and low BH masses, the scaling index a between the residence 
time and the period t^es oc f^^.^ will probe the physics of the circumbinary accretion disk 
and viscous orbital decay. At shorter periods and higher masses (roughly, at torb < few 
weeks for M > 10^ Mq), the identification of a /var oc t^l^ power-law would confirm that 
the orbital decay is driven by GWs. This would amount to an indirect, statistical detection 
of GW-driven SMBHBs, independent of any direct detection of GWs by LISA. This would 
also confirm that circumbinary gas is present at small orbital radii and is being perturbed 
by the BHs - and would thus serve as a proof of concept for finding LISA electromagnetic 
counterparts. 

The rest of this paper is organized as follows. In § [21 we discuss the evolution of 
binaries with different masses and mass-ratios, embedded in a circumbinary gas disk. We 
describe simplified models for the disk and for the binary-disk interaction, and emphasize 
that the binaries probe the distinct physical regimes in the disk, before GWs take over 
their evolution. In § [3l we discuss the possibility of searching for a population of coalescing 
SMBHBs among a catalog of luminous quasars, either based on their variability, or on shifts 
of their spectral lines. We discuss modest constraints available from existing surveys, and 
comment on specific recently detected individual SMBHB candidates. We then quantify the 
requirements for a detection in a future survey. In § |U we briefly summarize our results and 
offer our conclusions. When necessary in this paper, we adopt the background cos mological 



parameters Qm = 0.3, Q/y = 0.7, and Hq = 70 km s ^ Mpc ^ (IDunkley et al.ll2009l ). 



•^Throughout this manuscript, we will use the term "viscosity-driven evolution" to refer to the exchange 
of angular momentum and energy in the binary-disk system that arises from the combination of gas viscosity 
and the tidal torques from the binary. 
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2. Binary Evolution 



In this section, we describe the evolution of the orbital separation of a SMBH bi- 
nary. The basic picture we adopt is that the binary is embed ded in a thin circumbinary 



disk, with the plane of the disk aligii ed with the binary's orbit (IBardeen fc Peterson! Il975 



Ivanov. Papaloizou. &: Polnarevlll999l ). Initially, the orbital decay is dominated by viscous 
angular momentum exchange with the gas disk. However, the time-scale for viscous de- 
cay decreases relatively slowly as the orbital separation R decreases (tres oc 

^1/2 _ ^11/4 

see below) whereas the time-scale to decay due to gravitational radiation decreases steeply 
(tow oc R^). Therefore, generically, there exists a critical orbital radius -Rcrit, below which 
the decay is dominated by gravitational radiation. 

To describe the evolution quantitatively, we make several simplifying assumptions. 
The circumbinary gas is assumed to form a stan dard geometrically thin, o ptically thick, 
radiatively efficient, steady-state accretion disk ( Shakura fc SunyaevI 1973 ). We a ssume 



zero ecc entricity for both the binary and for the disk (justified by 



ever see 



Armitage fc Nataraian 



Hayasaki et al.l 



2007; 



2005 



Cuadra et al 



MacFadven k Milosavlievid l2008l : iDotti et al.l 12008 



Dotti et al 



200 



^, how- 



20091). and we assume co-planarity betwee n the disk and 



the binary ( Bardeen fc Peterson 19751 : Ivanov. Papaloizou. &: PolnarevI Il999l ). All of these 



assumptions may fail in the late stages of the merger (even before GW-driven decay be- 
gins). However, under these assumptions, the disk structure and the orbital decay have 
simple limiting power-law solutions, with the power-law indices depending on the choice for 
the underlying physics. These solutions are useful to describe the possible evolution of the 
binary, and to illustrate the point that the decay rate is generically a different - much flatter 
- function of torb than the tew ^ ^orb behavior in the GW-driven case. 

We emphasize that our aim here is not to provide accurate, self-consistent solutions 
for the co-evolution of the SMBH binary and circumbinary disk. Rather, we derive only 
gross scaling laws in various regimes - our main point is that these regimes and associated 
uncertainties, which are large, can in principle be probed observationally. 



2.1. Notation 



We ad opt the follow i ng no tation throughout this paper. We refer the reader to lShapiro fc Teukolsky 



(119861 ) and iFrank et al.l (120021 ) for general introductions to accretion disks. 



Physical constants: G is the gravitational constant; c is the speed of light; ke is 
the Boltzmann constant; ctsb is the Stefan-Boltzmann constant; ctt is the Thompson 
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cross section; /ig = riem-a/p is the mean mass per electron in units of hydrogen atom 
mass, mn, which satisfies /le = (1 + -^h)/2 for a fully ionized gas of both hydrogen 
and helium; is the mass fraction of hydrogen; /xq = 2/(3Xh + 1) is the mean 
molecular weight; k^s = /^eO"T/''^H is the electron scattering opacity; and Kfj = (8 x 
10^^ cm^ g" -M/iJp /(g cm-^lKT/K)-^/^ is the Rosseland mean absorption opacity in the 



free- free regime (jPadmanabhanll2002l : iRybicki fc Lightmanlll986l : IShapiro fc Teukolsky 



19861). 



BH parameters: Mi and M2 are the individual BH masses; M = Mi + M2 is the 
total BH mass; q = M2/M1 < 1 is the mass ratio; Qg = 4g/(l + g)^ is the normalized 
symmetric mass ratio; n = QsM/A is the reduced mass; R is the binary separation; 
Rq = R/ (1 + q) is the location of the lower-mass secondary, measured from the center 
of mass of the binary; Rs = 2GM/c^ is the Schwarzschild radius corresponding to 
the total mass; LEdd = 4:7[Gck~J^M is the Eddington luminosity for a BH of mass M; 
^Edd = -^Edd/(ec^) IS the Eddiugtou accretion rate with a radiative efficiency e; and 
^Edd = M/ Msdd = i^es^c/ (47rG) = (3.94 x 10'' yr) x fie^o.i is the characteristic time-scale 
associated with Eddington accretion. 

Disk parameters: H is the vertical scale height (the effective geometrical semi-thickness 
of the disk); p is the volumic gas density; S = p/{2H) is the surface density; Pgas is the 
gas pressure; Prad is the radiation pressure; P = Prad + Pgas is the total pressure; /? = 
-Pgas /(Prad + Pgas); T is the (midplauc) gas temperature; Teg is the effective temperature 
defined such that the locally emitted flux through an infinitesimal disk surface element 
is asBT^s', ^{r) is the Keplerian orbital angular velocity; Pa is the outer radius of 
the gap in the punctured circumbinary disk, measured from the center of mass of the 
binary; 77 is the anomalous dynamical viscosity; u = rj/p is the anomalous kinematic 
viscosity; a is the standard viscosity parameter of thin accretion disks; 6 is a constant, 
either or 1, determining whether viscosity scales with the total or just the gas pressure, 
so that rj = aP(3^VL~^\ k is the opacity of the disk material; r = (1/2)kS is the vertical 
optical depth; and fx is a constant defined such that fx = t'^T'^/T^q. Quantities 
with a A subscript (e.g. Qx, Sa, Hx) denote parameters in a steady-state disk around 
a single unperturbed accreting BH, computed at the radius Pa- Similarly, quantities 
with a subscript are evaluated at the position of the secondary Rq. 



With the above definitions, we proceed to define the dimensionless quantities r = R/Rs, 
rg = r/10^ M7 = M/(1O^M0), rh = M/MEdd, mo.i = m/0.1, ao.3 = a/0.3, eo.i = e/0.1, A = 
Rx/R, Kes = i^/kss, and kn = k/kq. Note that radii are measured from the center of mass of 
the binary. We adopt the set of fiducial values Xh = 0.75, pe = 0.875, po = 0.615, Jt = 3/4, 
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and denote values relative to these fiducial choices with a tilde, e.g. jle = /ie/0.875. Our 
fiducial binary+disk model is therefore chosen to he fie = P^o = fx = Q = Qs = Mj = mo.i = 
eo.i = ao.s = A = 1. Since all of our expressions can be written as products of power-laws in 
the physical parameters, the resulting expressions become tractable in these units. 



2.2. Thin Disk Models 



We next collect the basic expressions from the literature for accretion disk models under 
different physical conditions. We quot e the equations for a rang e of different steady thin 
disks, valid for a single accreting BH (IShakura &: SunyaevI Il973l ). We distinguish several 
cases: (i) whether the radiation or gas pressure provides the dominant vertical support, (ii) 
whether the opacity is dominated by electron scattering, Hes, or free-free absorption, Kg, and 
(Hi) whether the viscosity rj is proportional to the total pressure or the gas pressure (also 
known as a and P disk models, respectively). Based on these choices, the accretion disk can 
be divided radially into three distinct regions (IShapiro fc TeukolskyI 119861 ): 



1. Inner region: Radiation pressure and electron-scattering opacity dominate, P ~ -Prad, 
Kes ~ 1, valid inside -C rf^^^^'^ where rf^^^^'^ is defined in equations ffT2]) and ffT51) 
below. 

2. Middle region: Gas pressure and electron-scattering opacity dominate, P ~ Pgas, i^cs ~ 
1, valid between rf^^^^^ <^ < rg^^^, where r^^^^ is defined in equation f|T4l) below. 

3. Outer region: Gas pressure and free-free opacity dominate, P ~ -Pgas, % ~ 1, valid 
outside of rs >,rl^^^. 



In region (1), it makes a difference whether the viscosity is proportional to the total 
pressure or just the gas pressure, labeled below by 6 = or 1, (i.e. a or /5 disk) respectively. 
In all cas es, we assume tha t th e disk is optically thick, i .e. r ^ 1. We obtain S(r) and H{r) 
following ICoodmanJ fl2003h or ICoodman fc Tan I J2004h . 



E(r) = f^y^^^a-^/'n-y^M'/'^^^^ 



"^The choice fx = 3/4 is appropriate for a one -zone model where all the energy is dissipated near the 
midp lane and the opacity is con stant yertic a lly ( Shapiro fc Teukolsky I98(i Frank et al. 20021 : Armitagj 
2007 ). For reference, we note that Goodman ( 20031 ) and Sirko fc GoodmanI (|2003l) adopt different yalues of 
/t = 1 and /t = 3/8, respectively. 
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where 6 = or 1, and the radial dependence is imphcit in VL and [3. Here, /3(r) = Pgas/ (-frad + 
Pgas) which satisfies 

/3(l/2)+(l/10)(fe-l) / K \ 2/5 

1-/3 V/^omH; ^ ^ 

The asymptotic hmits of equations ^ and ([2]) can be obtained in regions (1-3), using 
equation ([3]). The resuhs are 

Inner region: 

/ • \ 3/5 

E(r) = (1.63xl0^gcm->f/.;^/^^;VVT-^a^^^^^ M^/^^/^f b = 1, (4) 

-,4„^^-2n, -l~-2 Z' ^ ^3/2 



= (2.50xl0^gcm-^)/i;iK-VT «a3(^— J ^a' if b = 0, (5) 

Tfl 

H{r) = {10.0Rs)fT — for arbitrary b. (6) 

eo.i 

Middle region: 

/ ■ \ 3/5 

S(r) = (1.63xl0^gcm-);.f;.-/^^-VV-a-f M^/V3-^/^ (7) 

Hir) = i3AlRs)^.:^/^W^'kU^'fTao^^^^^ (8) 
Outer region: 

E(r) = (2.61xl0^gcm-);.-/Vf%^/^V.-"^/"«-f (^) M,^/V3-/^ (9) 

/ . N 3/20 

Hir) = (3.08i?,)^-/io^-3/8~i/20^i43/i60^-Vio^^j M-/^Vf. (10) 

The boundaries between the inner/middle and middle/outer regions can be found from 
equations ([I])-(l3]), by requiring Pgas = -Prad and Kg = Kes, respectively. Note that Kff(r) oc 
pT'^l'^ depends on radius implicitl y through the density a nd the temperature. Using the 



(mid-plane) temperature given by iGoodman fc Tan I (12004 ). 



T(r) = (16772)-^/^^) a-i/5«:i/^M2/5fi3/5^-(V5)('>-i)^ (H) 
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we find that the transitions are located at the radii 

^gas/rad ^ Q AS2 kli' al{^\m,,, / e,.^^^^^^^ ifb = l, (12) 

= 0.515 /ir/ie/'''^L/'«ofK.i/eo.i)^'/^^M?^^^ if b = 0, (13) 



r 



■-/^ - 4.10/i-^/7T'^''(^ff/^es)-^/^(mo.i/eo.i)^/^ (14) 



3 



Note that the middle and outer regions differ only in their opacity laws, and the equa- 
tions in these two regions are equivalent (this can be seen by setting k^s = KfjKfj(r)/Kcs)- 
Since S, H, p, and T scale with a low power of kg, the radial dependence ends up being sim- 
ilar in the middle and outer regions. The distinction between these equations is nevertheless 
useful, since we can assume that kes — > 1 and kg —>■ 1 are constants in the middle and outer 
regions, respectively. 

We emphasize that equations (HI)- (fTOl) represent only a very non-exhaustive subset of 
solutions even for radiatively efficient steady thin accretion disks. In particular, at large 
radii, there are several effects that can invalidate the disk model described by these equations. 
First, these solutions assume that the self-gravity of the disk is negligible. This assumption 
becomes invalid at radii where the Toomre Q-parameter equals unity, 

rf = 12.6/i-^/^/i-/^^^°/^^L/^aJ4^(mo.^^ if /^e. - 1 (15) 

rf = 30.99/io-^/if f'^''/''kl/''alT (^o.i/eo.i)"''/'' M^'^'' if k^ 1. (16) 

Beyond these radii, the disk is commonly believed to be un stable to fragmentation . 



Second, at large radii, the disk can also become optically thin (see lSirko fc Goodman! 12003 
where solutions can be obtained by fixing the Toomre parameter in the outermost region 
at Q = 1). At th e se binary separations, the disks ma y not actually be geometrically thin 



( jPotti et al.l I2OO8I : iMayer. Kazantzidis. fc Escalal |2008| ) , and slim or thick solutions might 



instead be relevant. Third, beyond the radii where the disk temperature falls below ~ 
lO^K, the gas becomes neutral. The corresponding change in opacity will modify the disk 
structure, and the disk may become susceptible to ionization instabilities (although see 



Menou fc Quataertll200ll ). Finally, at large radii (where the orbital velocity > 100 km/s), 
the gravitational potential of the galaxy can no longer be ignored. These regimes, however, 
turn out to correspond to separations larger than we are interested in the present paper, for 
BH masses above ~ 10^ Mq (as will be shown in Figures [Hand [2] below). 



2.2.1. Comparison with Other Results 



We have verified our solutions numerically by su bstituting them back i nto the fundamen- 
tal conservation equations of thin accretion disks (jshapiro fc Teukolsky 1986 : Frank et al, 
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20021 ). Moreover, equations dl])-® agree with those quoted in iGoodman fc Tan 



Equations (I5l)-(l6]) are also consistent with lShapiro fc TeukolskyI (119861 . page 441), for the 
P = Prad, 1^ = /Ces, 6 = model. It is also reassuring that equations (l^- (fTU|) are consistent 
with those in Frank et al. (2002, Sec. 8.1, p. 244)]^ Also note that, owing to the wea k 
dependence on % , our numerical factors are very similar to those in iFrank et al.l (120021 ) , 
even though is defined to be two orders of magnitude larger there than the value we 
adopted here (to be consistent with most other textbooks). 



2.3. Binary — Disk Evolution 

Here we collect and summarize the most important formulae describing the interaction 
between a binary and the accretion disk in order to identify the mechanism that drives the 
orbital decay of the binary during the final stages of the merger, as a function of binary 
separation (the choices being GW driven inspiral and tidal- viscous torques). This will allow 
us to explicitly compute the residence time tres = —R{dR/dt)~^ that an individual binary 
spends at each orbital separation R, or at the corresponding orbital timescale torb- 

The formulae collected in this section will also allow us to quantify the binary separation 
at which the viscous evolution of the disk is decoupled from the increasingly rapid, GW- 
driven orbital decay of the binary. We provide results for a and /?-disks, and give analytic 
results function of binary and disk parameters. 



2.3.1. Disk- versus Secondary-Dominated Orbital Decay 



In general, the evolution of a SMBH binary in a thin disk is analogous to planetary 
migration (see, e.g. lArmitagdl2007l ) . In the limit of a very low-mass companion (g -C 1), the 
interaction between the planet and the disk is linear. In addition to co-rotation resonances, 
the density waves excited in the gas at discrete Lindblad resonances with the binary exert a 
large net torque on the binary, leading to rapid, so-called Type-I migration, which occurs on 



^It appears that I Goodman &: Tan I (|2004l ) contain the foUowing typographical errors: p oc 
Q,4/5^-4/5^j4/5^-i8/5 ^^^^^ should be p cx a'^^'^ M^^^^r^^^^^ , and c, cx K^^'^lEddM^r~^^'^° 

in their eq. 21 should be Cg 
regions, respectively. 



Kl-EddM^r^ ^^"^ and Cs oc K^/^^/g'^jAf-, ^^^^r^ ^^"^^ in the inner and middle 



^However, there appears to be a typographical error in their quoted scaling H oc , which should 

- 2 /2O 

instead read as iJ oc Afjg i so that S = 2pH is satisfied for all M26. 
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a time-scale much sh orter than the local viscous time-scale (e.g.. lTanaka. Takeuchi. fc Ward 
2002l : lArmitag3l2007h . 



If the binary is massive enough for the tidal torque to dominate over the viscous torque 
in the disk, the interaction becomes non-linear, and a gap is opened in the disk, extending 
to the outer radius R\ = \R. The condition for a gap to open is that the mass ratio exceeds 



the critical value q > max{(i7o/^o) , {lOay'^Ho/ Rof'^} (e.g. lR,afikovl 120021 . note that Hq 



is evaluated at the position of the secondary Rq). For binaries that are not in the GW- 
driven regime, and for which the disk mass exceeds the mass of the secondary (see below), 
this typically translates into the very modest requirement q > 10~^. This is satisfied for 
all SMBH binaries that may produce the electromagnetic signatures we discuss below. The 
exceptions are the so-called extreme mass-ratio binary inspirals (EMRI's) with q < 10^"^ 
(i.e. a stellar-mass object coalescing with a SMBH). In this paper, we focus on SMBH 
binaries, and therefore in the rest of this paper, we neglect Type I migration. 

If the secondary's mass satisfies the above gap-opening threshold, but is still small 
compared to the local disk mass, then it acts as an angular momentum bridge for the disk, 
and the secondary's orbital evolution is simply determined by the viscous diffusion time. 



tu 



r 
r,. 



3 yo 



2n- 



0^0 



M 



(17) 



where we have used M = 37n/nSn which follow s directly from angular momentum conserva- 
tion in steady disks (IShakura fc Sunyaevlll973l ). The orbital decay of the binary in this limit 
is analogous to disk-dominated Type-II planetary migration. 

In practice, the assumption that the local disk mass exceeds the secondary's mass often 
fails. In this case, analogous to "planet-dominated" Type-II migration, the angular mo- 
mentum of the binary can still be absorbed by the gaseous disk outside the gap, and the 
viscosity of the gas can drive the binary toward merger. However, migration is slower, and 
the time-scale in this regime, t^, is longer than t^. An estimate of the slowing factor is g^*', 
where 



4:71 RiE 



Qb 



0^0 



/i 



2M 



-t„ 



8m ty 

Qs ^Edd 



is a measure of the lack of local disk-mass dominance (ISyer &: Clarkd Il995l . but note that 
our is denoted by "i?" in their original definition), which is less than unity in this case, 
and A; is a constant defined as 



1-1 



81nS 
d\nM 



-1 









if qe < 1 & Kes ^ 1 
if Qb < 1 & ^fT ^ 1 
if qB > 1 



(19) 
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Thus, the separation of the binary in this case is driven inward on the timescale 

ts = -- = qB-'U=(^^-^Xu ifqB<l. (20) 



Note tliat the viscous time-scale ti, in disk-dominated hmit (eq. [T71) should be evaluated 
at the position of the secondary Rq = R/ (1 + q), while the quantities entering the time-scale 
ts for the secondary-dominated type-II migration of more mass ive binaries ((/r < 1. eg. [201). 



should be evaluated at the outer edge of the cavity, Rx = XR (IMacFadyen &: Milosavljevic 



20081 ). In order to avoid a discontinuous jump in the migration time-scale at the qb = I 



transition, below we will omit this distinction, and evaluate both time-scales at Rx = XR. 



2.3.2. Evolution of Individual Binaries From Large to Small Radii 

Using the steady thin disk model outlined above, we can calculate the rate at which the 
binary is driven inward by the gas. We will also estimate the rate at which the inner edge 
of the punctured gaseous disk follows the binary due to its viscosity. From the preceding 
discussion, we see that both the viscous time-scale and the orbital decay rate depend on 
whether the binary is located in the inner /middle/outer region of the disk; and also on 
whether the local disk mass is larger /smaller than the mass of the smaller SMBH. For 
completeness, we here obtain and quote the residence time as a function of orbital radius 
and orbital time, in each of these 3x2 = 6 regimes. We then construct the self-consistent 
evolution of individual binaries, with different masses and mass-ratios, across the relevant 
regimes. 

We first consider the timescale t^, and assume that the secondary perturbs the disk at 
the radius Rx = XR. This is the relevant regime initially, at large binary separations, when 
the disk mass enclosed within the secondary's orbit is large. In this regime, we find, 

t, = (2.82 X lOVr) X K,-7f ^ao:^ (^^^ \rX'/^/' if b = and rg < rr/''"^(21) 



= (5.96 X 10^ yr) x f^l^' f^t^'k-'^' f^WJ' (^) M?/^A^/Vj/^ 

ifrr/-'^«r3<rr/", 
or if b = 1 and < r^'^'''''^'^'^ 

/ . \ -3/10 

(7.37 X lOVr) X /i^^/if ^/^7t '''/'°«af ( ^ ) M?/^A^/Vf 



(22) 



, eo.i 

cs/ff 



ifr3>rr^ (23) 
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eo.i / J 

(26) 

-3/10 / , N 5/6 



The above can be expressed as a function of the orbital time of the binary, 

t^^^ = ^ = v^vr^r^/^ = 2.81 x lO'^rl/' = (0.88 M,rl^') yr, (24) 
which results in 

= (7.48 X 10^ yr) x k^^f^'a.l ' M^^'X'I^ 

if b = and rs < iT''""^ (25) 

/ • \ -2/5 / , \ 14/15 

U = (6.73 X lOVr) X ^y5^f ^-V^-^aoT f ^ M.^/^V/M ^ ) 

;f „gas/rad ^ „ ^ „cs/ff 

II rg ^ ^3 ^ rg , 
or if b = 1 and r3 < r|^*'^'^'**^ 

U = (8.21 X lOVr) X /i^^/if %'/^7T'''^'°«af " M.^^/^^A^/^ (^^) 

ifr3>rf^ (27) 

The measure of disk dominance can be calculated by substituting the viscous time-scale 
into equation f|T8l) . 

if b = 0andr3<rr/'^'''' (28) 

/ • \ 3/5 

= 0.011 ;i-/^/if^-/5/-2^-4/5^^^ M,^/^-A^/Vl/^ 

ifrr/^^'«r3<rf^^^^^^ ^29) 
or if b = 1 and r3 < r|'^'''^'^'^'^ 

= (1.49 X 10-) ^-/^/if %^/^7.-"^/"«o1/^ ) M,^/^-A^/V3^/^ 

ifr3>rf^. (30) 

In order to decide whether the evolution indeed follows the "disk-dominated" decay (on the 
viscous timescale ty) or the secondary-dominated decay (on the longer time-scale tg), one 
should examine whether > 1 or < 1 is satisfied, respectively. From equations (!28|) - (!30|) . 
we find that the transition occurs at 



/ \ 2/7 

- 3.61/ire/^/WJr^ M-V^A-^ 

eo.i 
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if b = 0andr3<rr/'""'' (31) 

/ • \ -3/7 



■ r gas/rad , es/ff 

if rf ' < 13 < r ' 
or if b = 1 and 13 < ig 



gas / rad 

(32) 



/ . N -14/25 

3^/^ = 182/i^^/2^/i-'/^4/^7T''^'°Var(^) Mf'!^\f^\-^ 



, eo.i 

es/ff 



ifr3>rr", (33) 

Note that with the exception of very unequal masses q < 0.01, the transition takes place 
well in the outer region of the disk, with rij^ > 10. At smaller radii, the binary is driven 
viscously on the timescale (rather than t,^). 

The "secondary — dominated" Type-II decay timescales relevant at these radii can be 
obtained by substituting equations (!2T!) - (l23l) into equation ( !20!l 

/ . N -13/8 

4 = (6.15xl0Vr)x/i3/8^-//^/^-^/X:f M',l\ll^\^^l'\f^' 



. eo.i , 

if b = andrg <rr/''"'' (34) 

/ • \ -5/8 

4 = (7.40xl0^yr)x^V2^y--V8/-5/4^-V^ M,^/^g3/8^r/8,r/8 



\t „gas/rad ^ „ ^ „es/ff 

It rg <C rs ^ , 
or if b = 1 and rs < ^^^^'^'^ 



(35) 



/ . N -10/17 

ts = (1.07 X lOVr) X /i^A7/i^^/34^-Vi7-i43/i36^-8/n Mf^J/i^A^^/^Vf/^^^ 

ifr3>rr^^ (36) 
or, in terms of torb using equation (|2^ . 

/ . \ -13/8 / , \ 35/24 

4 = (3.60xl0Vr)x/i3/8^-5/^/-V4^-f Mf^VA^'/'M 

if b = 0andr3<rr/''"'' (37) 

/ • \ -5/8 / , N 7/12 

= (7.98xl0^yr)x;iV2-/--V8/-^^^^^ M^^gf^AV^ j 



.f gas/rad es/ff 

if r| ' < rg < rg ' , 

or if b = 1 and rg < ^^^^'^'^ 



(38) 
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^6,_^ -9/17 ~15/34~ -1/17 ?-143/136^ 



;i.l4 X lO^r) X /i^/"/io 

.f ^ es/flf 

II rs ^ fs ' 



8/17 
0.3 



^0.1 

eo.i 



-10/17 



-11/51 7/17^25/34 / torh \ 



25/51 



(39) 



Finally, at a still smaller radius, the orbital decay will be dominated by gravitational 
wave emission. The GW-driven decay timescale in the leading order (Newtonian) approxi- 
mation, is 



t 



GW 



5Rs _ 



V4 



:i.ll X lOVr) X g;iM-'/' (^-^^ 



8/3 



(40) 



This approximation is adequate for our purposes, since post-Newtonian corre ctions do not 
become appreciable until the final ~ day of the merger (see, e.g.. Figure 5 in iKocsis et al. 



20081 ) . Note that tow defined above differs from the total time to merger, (defined as the 
binary separation decreasing to zero), which is often used in the literature, and which occurs 
at t^^'^^ = tQ-^/i. What is the radius at which tew becomes smaller than the time-scale for 
Type-II orbital decay? Let us express this transition in terms of the radius that satisfies 
^GW = Pcw/sts, where Pgw/s is a fixed constant of order unity: 



s/GW 



„s/GW 



20/29 -10/29 
'^0.3 



0.587/i^/29^;^20/29j- 

if b = and < r|'^'^'''^'^'^ 



m 
eo.i 



-26/29 



^-6/29^22/29;^35/29^16/29 



GW/s 



(41) 



0.470 /i^/-/2S/^^^-/-/,-^/^a-f ^ ( ^ ) M-/-g-/25^7/25^8/^5^^ 



eo.i 



-1/5 



if r' 



gas / rad 



< r3 ^ r3 



es/ff 



s/GW 



or if b = 1 and < x^^^^^'^ 



0.545 /i^/37/if'^^2/"7T'''^'''-"''^'" 



(42) 



a, 



0.3 



eo.i 



-20/111 



11 
GW/s 



• r ^ cs/ff 
if ^3 ^ fs • 



^-10/111^16/37^25/111^34A 

(43) 



The corresponding critical radius is around ~ 500i?5 for system parameters near the assumed 
fiducial values. The critical radius, however, is significantly closer in for very massive, and 
very unequal-mass binaries (i.e. for M = 10^ Mq and q = 0.01; see Fig. H] below). Interest- 
ingly, the critical radius is quite insensitive to the BH mass and accretion rate (i.e. to Mj 
and m). Note that the viscous timescale, t^, describing gas accretion, is faster than tg, which 
indicates that at the time when GW starts driving the evolution, the viscous inward diffusion 
of gas can initially still follow the binary. However, the comparison of equations fl2T|) - fl23|) and 
equation fHOj) shows that as the binary orbit shrinks further, the viscous time-scale always 
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decreases less rapidly than the GW inspiral timescale, so that eventually the evolution of 
the gaseous disk will decouple from that of the binary. Let us find the critical radius, r'^^'^^ ^ 
where GW inspiral outpaces viscous gas accretion. We find that in most cases, this critical 
radius is not relevant for the orbital decay of the BHs themselves, because the transition to 
secondary-driven orbital decay always takes place before GWs start dominating the decayJil 
However, this critical radius is relevant for the behavior of the disk: it provides an estimate 
for the time when the punctured disk decouples from the GW-driven binary, and effectively 
stops evolving (and also for the size of the inner gap at this time and onward). By requiring 
^GW = Pcw/utu, where Pgw/u is a constant coefficient of order unity, we obtain: 



0.202 ^-^^ 



0.3 



m 



eo.i 



GW/!/ 



if b 



0.222 /iVi3-/i3~- 



1/13 /-10/13 



4/13 
0.3 



• r gas/rad , es/ff 

if r^ ' < rg < rg ' , 

or if b = 1 and rg < x^^^^^ 



0.183 /iy^^/if^^-/^V. 

• r ^ cs/ff 

if ^3 ^ ' . 



Ot, 



0.3 



^0.1 

eo.i 



13/20 -16/55 / "^0.1 



eo.i 



and rs < x^^^^^ 



-2/13 



^1/13^5/13^/13^5/^3^^ 



-6/55 



^4/55-4/11^5/11^4/n 



(44) 



(45) 



GW/i/ 



(46) 



The appropriate choices for /3gw/s and Pgw/u are poorly known, but /?gw/s may be 
reasonably taken to be ~ 1 when the binary is first driven by GW emission, rather than 
by tidal interaction with the ga s. The simplest choice for Pgw/,,, adopte d in many previous 



studies, is also Pgw/u = 1 (e.g. lArmitage fc Natarajanll2005l : iLoebl 120071 ). However, the gas 



infiow rate across the edge of the central gap will be increased due to the steep density 
and press ure gradient (ILvnden-Bell fc Pringldll974j ). which will delay the decoupling. This 
motivated iMilosavljevic fc Phinney J2005h to adopt Pgw/u ~ 0.1 (the value describing the 
limiting case of an infinitely sharp edge). 

Adopting Pgw/s = 1 in equations (HTi)-(H3l) then yields the radius where the binary evolu- 
tion changes from being viscosity-driven to GW-driven, and Pgw/u = 0.1 in equations flUl) - 
(146|) gives the separation at which the disk totally decouples from t he binary and the radius 
of the gap "freezes" . These expressions generalize the results of IMilosavljevic fc Phinney 
(120051), who restricted their analysis to the 6 = 1 case, and focused on the behavior of 



^The exceptions to this are the most-massive, M > IQi" M0, equal-mass binaries, and only if 6 = 1 
is assumed - in this case, the GW-inspiral takes over in a radiation-pressure dominated disk, in the disk- 
dominated regime, i.e. before the transition to the secondary-dominated regime. 
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the disk at decoupling, rather than the orbital evolution of the binary. In particular, 



Milosavljevic fc Phinneyl (120051 ) evaluate disk conditions at the single radius at the edge 
of the gap, at the time of decoupling, and do not discuss the transition from the disk- to 
the secondary-dominated decay, or other details of the binary's orbital decay. The binary 
separation at decoupling is of order 7-^"^'^/^^ ^ Q.l for both the gas pressure dominated 
models and the radiation pressure dominated case with 6 = 1. In these cases, the transition 
between viscosity and GW-driven decay and the disk decoupling take place in relatively 
quick succession, since 7-3'^'^^'"^^ depends weakly on (3gw/u- The delay between these two 
events is much longer for the radiation pressure dominated regime when 6 = 0, since in this 
case the viscosity, which is proportional to the total, rather than just the gas pressure, is 
much larger, and the gas can follow the binary nearly all the way to merger (at least for 
large m). In this case, the result is also extremely sensitive to the accretion rate and the 
binary mass ratio. Generically, for a fixed total binary mass, the decoupling occurs at the 
largest separations for nearly equal masses. 

Interestingly, the decay rate of a given individual binary can decelerate and accelerate 
during its evolution, according to the variations in the local disk environment at each instan- 
taneous binary separation. The evolutionary tracks of binaries with four different choices 
for the total mass (M = 10'^, 10^, 10^, and 10^ Mq) and two different mass ratios (g = 1 
and q = 0.01) are shown in Figures [1] and [2l In both figures, we assume that the viscosity 
is proportional to the total pressure (6 = 0). The motivation for this choice is to illustrate 
the effect of the additional radiation pressure-related viscosity on the orbital decay (which 
is not present in the 6 = 1 case). We note t hat a phenomenological 6 = disk is known 



to suffer from a formal thermal instability (e.g. iLightman fc Eardleylll974l ): recent magneto 



hydrodynamical simulations, however, found such disks ther mally stable (while accountin g 



for the correlation between viscosity and radiation pressure iHirose. Krolik fc Blaed l2009l ). 
These figures show the residence time as a function of the orbital time. They demonstrate 
that the evolution of the binary in most cases proceeds through the following distinct stages. 

(i) Disk-dominated viscous evolution. Initially, at large separations (shown in blue 
curves), the binary is strongly coupled to the circumbinary disk and evolves on the viscous 
time-scale tvisc (analogous to "disk-dominated" planetary migration). The radius of the 
gap follows the binary. During this stage, tres ~ is proportional to r^l^-r^l'^ (the range 
corresponding to the choice 6 = vs. 6 = 1) for radiation pressure, or 7-''/5_^5/4 i^t^,^ 
range corresponding to the choice of dominant opacity being electron scattering or free-free 
absorption) for gas pressure dominated disks (see eqs. [2IH2SD- These decay rates translate 
into ty oc if^orb^^^^rb' ^orb^^^^orb' cascs respectively (see eqs. [^5H27|) . Note, 

however, that for nearly equal-mass binaries, the separations have to be quite large to 
correspond to this disk-dominated regime - falling into the outer regions of the disk, which 
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are unstable to fragmentation (the orbital radii where the disks are marginally Toomre- 
stable are marked with large dots). Therefore, depending on the behavior of the gas disk 
beyond this radius, this early stage of disk-dominated viscous evolution may exist only for 
unequal-mass binaries. As shown in Figure [H disk-dominated viscous evolution may be 
realized in a stable disk for binaries with M < 10^ Mq and q ~ 0.01; in these cases, the 
binaries are in the free-free opacity and gas-pressure dominated regions of the disk, so the 
relevant scaling is tj, oc r^^* oc t^J^^. 

(a) Secondary-dominated viscous evolution. As the binary separation shrinks below 
R^h ^ i{)^Rg {l^^Rs) for mass ratios g ~ 1 (g ~ 0.01), the binary mass starts to dominate 
over the local disk mass, and the binary evolves more slowly, according to "secondary- 
dominated" decay (analogous to "planet-dominated" Type-II migration). During this stage, 
the GW emission is still negligible, and the decay time-scale can be obtained from equa- 
tions (!MI) - (!36l) . and t^ oc r'^l^-r^^/^^ for radiation pressure (with 6 = — 1), or 7-'i'/8_j,25/34 
gas pressure dominated (with electron scattering vs. free-free opacity) disks, implying that 
^ ^oJo^^^orh^i ^orb^^^orb^^ ' cascs respectively (see eqs. |37H39|) . As can be 

seen from Figures [1] and Ej on orbital time-scales between weeks to years, each of these scal- 
ings is relevant for some choice of binary masses. However, the transition to GW-domination 
always takes place either in the "inner" or "middle" disk region. 

(Hi) GW-dominated evolution. Still later, within the radius R^l*^^ ~ 500i?5 for systems 
with parameters close to the fiducial values, the binary's orbital evolution starts to be driven 
primarily by GWs, but the outer edge of the gap can still diffuse inward and follow the 
binary. During this stage, the decay time-scale is tow oc oc t^^j^. 

(iv) Gas disk decoupled. Finally, within R'g'^^ ~ IOO-R5 the binary is entirely driven 
by GWs and the binary falls in much more quickly than the outer edge of the gap is able to 
move inward. 

The above ordering of events is valid for a broad range of binary and disk parameters. 
Note that the ultimate fate of the gas inside the b inary's orbit is left unspecified in our 



considerations (see, e.g. lArmitage &: Natarajanll2002l . for a possible outcome) 



In addition to the above sequence of events describing the evolution of individual bina- 
ries, several interesting conclusions can be drawn from Figures [1] and [2l 

1. Goalescing binaries have a non-negligible abundance. First, binaries with masses in 
the range 10^ — IO^Mq may be both bright and common enough to be detectable in a 
survey, provided they have bright emission. Indeed, Figures [1] and [2] show that these 
binaries spend a non-negligible fraction (> 10^'^) of their total fiducial lifetime of 10^ 
years at orbital time-scales between 1 day < torb ^ 1 year (the total lifetime will be 
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justified below). It is feasible, in principle, to look for variability on these time-scales, 
and the residence times shown on the figures suggest that these variables may not be 
uncommon among bright AGN. We will discuss this possibility further in § 13.11 below. 

2. Disk- and GW-driven evolution may both be observationally relevant. Figures [1] and 
[2] also show that the transition from gas- to GW-driven evolution can occur within 
this "observational window" . For example, at the fixed orbital time-scale of torb = 10 
weeks, equal-mass binaries above lO^M© are GW-driven, and below this mass, they 
are gas-driven. 

3. Secondary-dominated evolution cannot be ignored. Essentially all binaries at the orbital 
times relevant for actual surveys (again, between 1 day < torb ^ 1 year) that are 
gas-driven are in the regime of "secondary-dominated" type-II orbital decay (referred 
to as stage (ii) above). Likewise, the transition from "gas-driven" to "GW-driven" 
evolution always occurs from the "secondary-dominated" type-II decay regime. In 
previous works wh ose primary focus was on the behavior of ga s at (and after) the time 



of decoupling (e.g. iMilosavljevic &: Phinneyll2005l : iLoebl 120071 ). this intermediary step 



which is important for the orbital decay of the binary, is not discussed. 

4. Observed binaries could probe all three disk regions. Interestingly, among the 10^ — 
lO^M© binaries with 1 day < torb ^ 1 year, it appears that all three of the disk 
regions (inner/middle/outer) enumerated in § 12.21 can be observationally relevant (i.e., 
gas-driven binaries can be found in each of these three disk regions). 

5. Viscous evolution is non-negligible even in the LISA regime. The comparison of Fig- 
ures [T] and [2] shows that unequal-mass binaries evolve more rapidly when they are 
gas-driven. Consequently, they make the transition to the GW-driven stage quite late 
in their evolution. In particular, binaries enter LISA^s detection range at the approx- 
imate observed GW frequency of few = 0.03 mHz. This corresponds to an observed 
orbital time (on Earth) of torb = 2//gw = 0.11 week. We find that at this orbital time, 
viscous evolution is not necessarily negligible. Figures [T] and [2] show that equal-mass 
binaries with M < 10^ Mq, and q = 0.01 binaries with M < 10^ are just making 
the transition to the GW-driven regime as they enter the LISA band. 

6. Total decay time in a stable disk is consistent with quasar lifetime. As Figures [T] and [2] 
show, the residence time at the radius at which Q = 1 is, in all cases, close to (although 
somewhat shorter) than the fiducial quasar lifetime of 10^ years. It is plausible that 
SMBHs become luminous, and act as quasars, only once they are embedded in stable 
circumbinary accretion disk. The fact that it takes ~ 10'' years for the binary to 
evolve from the outer edge of a stable disk to coalescence is therefore consistent with 
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the idea proposed in this paper, that there is a one-to-one correspondence between 
coalescing SMBHs and quasars (ahhough, as mentioned above, there are caveats that 
can invahdate the steady disk models at the relevant large radii). 



The possible implication of conclusion n o . 5 ab ove for LISA merits some further elab- 



oration. As discussed, e.g., in ISesana et al.l (120041 ). individual binaries can contribute to 



the LISA data stream in several ways. Sources can be divided into two types, based on 
whether they evolve significantly on a time-scale of ~ 3 years, the duration of the LISA 
experiment. Binaries caught at an orbital separation with short enough residence times for 
the frequency-evolution to be measurable are sometimes referred to as "gravitational sirens" 
or "gravitational inspirals" . Figures [1] and [2] show that during the last several years of the 
coalescence, the orbital evolution is always strongly GW-dominated, even for the lowest- 
mass BHs, and therefore the GW waveform of these rapidly evolving sources (including those 
whose actual coalescence is detected by LISA) will not be affected by the gas disk. 

Binaries that have a much longer residence time at some fixed frequency in LISA^s 
band represent "stationary" sources whose frequency remains roughly constant during the 
LISA mission lifetime. These sources could, in principle, be individually detectable by LISA. 
However, in practice, they are likely to accumulate sufficient signal-to-noise for detection 



only in the last few hundred years of their coalescence (see, e.g.. Figure 2 in ISesana et al. 



20051 for the detectability of g = 0.1 binaries as a function of their look-back time from the 
merger). Figures [T] and [2] show that viscous processes can significantly speed up the evolution 
of binaries only at somewhat larger look-back times (note that the look-back time is 4 times 
shorter than the residence time in the pure GW-driven case). The cumulative signal from a 
collection of faint stationary sources can, however, still add up to an unresolved background 
that is detectable, depending on the the cosmic evolution of the BH merger rate and the 
instrumental noise of LISA. The presence of the gas disks could reduce any such background 
that is present (compared to a prediction that assumes pure GW-driven evolution at LISA- 
frequencies) . 

In Figures [Hand [21 we have showed the evolution of the binary as a function of its orbital 
period. This will be particularly useful for assessing the detectability of such binaries in a 
survey for periodically variable sources (§ 13.11 below). In Figures [3] and [H we show, instead, 
the evolution of the same set of binaries, but as a function of their orbital separation. The x- 
axis on these figures is shown in units of Rs, with the corresponding orbital velocities shown 
by the labels on the top axis. This figure directly reveals that relatively more massive binaries 
(M > 10'' Mq) spend a significant time at orbital velocities of several thousand km s~^. Such 
orbital speeds may be detectable in the spectra of individual sources, providing an alternative 
to the detection based on periodic fiux variations (see § 13.41 below). 
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Finally, the conclusions enumerated above also highlight the large uncertainty in the res- 
idence times predicted in Figures [1] and [2|, caused by our idealized treatment of "secondary- 
dominated" type-II orbital decay. One immediate additional source of uncertainty is the 
choice of 6 = vs. 6=1. Before entering the GW-driven regime, most of the equal-mass 
binaries (Fig. [1]) are in the gas-pressure dominated region of the disk, but unequal-mass 
binaries (Fig. [2]) are in the radiation-pressure dominated region. Therefore, whether the 
viscosity is proportional to the total pressure or just the gas pressure makes little difference 
to the near-equal mass binaries. However, it makes a significant difference for unequal-mass 
binaries with M > 10^ M0. To show this explicitly, in Figure O the upper vs. lower curves 
contrast the evolution in the 6 = 1 vs 6 = case, respectively. As expected, once the binary 
approaches the radiation-pressure dominated regime, the evolution is significantly slower in 
the 6 = 1 case. The difference is most pronounced for the most massive (10^ M©) binary. 
For this system, the transition to GW-domination also occurs at a larger orbital time (~ 10^ 
weeks for 6 = 1, vs. 10 weeks for 6 = 0). 



2.4. Type II Orbital Decay in a Non-Steady Disk 

For simplicity, above we calculated the timescales and in steady thin disk models. 
However, as noted above, this highly idealize d mod el makes several crucial assumptions. 



In particular, llvanov. Papaloizou. &: PolnarevI (119991 . hereafter IPP) considered the tidal- 



viscous interaction of an unequal mass binary (g -C 1) with a time- dependent accretion disk. 
They assumed that the accretion disk is initially described by the steady-state solution for a 
single BH, and then considered the modifications due to tidal torques from a secondary BH. 
The torques are turned on suddenly at some moment to 5 when the secondary, whose mass is 
M2 = qMi ~ /i -C Ml, is at an orbital radius rg that encloses a disk mass M^q ^ ^2- The 
torques are assumed to be concentrated in a narrow ring near the secondary's orbit, which 
results in a pile-up of material near the outer edge of the disk cavity. They found (see their 
eq. 58) that this results in a decay time-scale of 



where M is the initial steady-state accretion rate, t^^ is the initial viscous time (at t = 
to), qBfl is the disk dominance parameter at t = to (see eq. [18]), = rt,(t) is the time- 
dependent position of the secondary and tq = rb(to) is its initial position, and rt/ro = 
[1 — ^S{t^^'^^''^^^'^'^^ — 1)]^, with the dimensionless time r = 2/3o{t/t,y), implying that 

^ 4c/(5c+b) 

1 +7*^ - Vrb/ro 



(4J 
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where S = {M / fi)ti^^o = qBfi/2, and a, b, and k are defined such that z/ = kT,°-r^, c = 
2(a + l)-6, = [c(2c + a)/[2(a+ l)]][(2c + a)/(2c+ 7 = 2c/[/5o(5c + 6)], 
and t,^ is the unperturbed viscous timescale given by equation (ITTI) and calculated exphcitly 
belowjfl If the opacity is dominated by electron scattering, then a = 2/3, 6 = 1, c = 7/3, 
j3o = 1.126, and 7 = 0.327, while for the free-free process a = 3/7, b = 15/14, c = 25/14, 
Po = 0.726, and 7 = 0.492. From equations (07]) and (08]) we find 

1/2 f >| 2{a+l)/{5c+fe) 

t:PP=f^V^) ^^^^1 ■ (49) 

\qB,oJ \roJ [2 + jqB -2y/rb/ro J 

Here, the radial-evolution given by equation is fat least initially ) not a simple power- 



law. Most importantly, as noted by llvanov. Papaloizou. fc PolnarevI (Il999l ). the pile-up of 



the disk material causes the binary decay to slow-down even more than estimated for a 
steady disk based on the "disk-dominance" parameter (eq. [20] above). 

In Figu res E] and [71 we illustrate the impact of allowing the disk to evolve. In the above 



approach of [ivanov. Papaloizou. fc Polnarev[ (Il999[ ). we have to specify when the interaction 



between the secondary and the disk is turned on. In Figure [6l we assume that the interaction 
begins as soon as the disk dominance parameter reaches qb = 1. In Figure [71 we delay the 
onset of the interaction to qb = 0.01. In both figures, the new dotted (magenta) curves 
denote the binary's residence time in the time-dependent disk. Note that in the latter case, 
in Figure [TJ the residence time undergoes a discrete jump when the disk-binary interaction 
is turned on: the binary stalls, and does not move initially, until the mass of material that 
has piled up is of the order of the secondary's mass. As the figures show, these residence 
times are indeed significantly longer than in the steady-disks. At relatively late times after 
the interaction is assumed to turn on, the residence times asymptote to power-law forms. 
Most significantly, for each of the binaries shown in Figures [6] and [TJ the transition to the 
GW-driven regime occurs significantly earlier due to the evolution of the disk; just before 
this transition to the GW-driven regime, the residence times are longer by ~ two orders of 
magnitude compared to a steady disk. 



3. Observational Search for SMBH Binaries Among Luminous AGN 

In the rest of this paper, we will discuss identifying coalescing SMBHBs with quasars, 
and interpreting the residence time tres as the duty-cycle for exhibiting periodic variability on 



^In IPP, refers to the standard gas pressure dominated accretion disk, which we generahze to radiation- 
dominated disks below. 
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the observed time-scale tvar ~ (1 + z)torh- Our broad justification for these hypotheses is the 
generic idea, advanced in numerou s other works, that quasars are activated in major galaxy 



Hopkins et al.ll2007al and references therein). S ince SMBHs are believed to be 



mergers (e.g. 

common in galactic nuclei (at least a t low redshifts z <3; see lMenou. Haiman fc Narayanan 



2OOII : iLippai. Frei fc Haimaru l2008bl and references therein) there could then arguably be a 



one-to-one correspondence between the quasar phenomenon and SMBHB coalescences. 



3.1. A Simplified Model for the Population of Periodic AGN 

The total quasar lifetime, defined as the cumulative duration (possibly over multiple 
episodes) for an individual source to produce bright emission near the Eddington limit, is 
generally belie ved to be tq ~ few x 10'^ yr, based on several lines of observational evidence 



(IMartinil I2OO4J ) . As discussed above, this value is consistent with the time-scale it takes for 
a binary SMBH to evolve to coalescence, starting from the outer edge of a gravitationally 
stable thin a-disk. Therefore, we hypothesize that the luminous quasar phase coincides with 
this last stage in the merger of the two SMBHs. Of course, it is possible that the quasar 
phase occurs either long before or after coalescence - in either case, there would be no bright 
emission to observe during the last stages, as hypothesized here@ 

We next assume that during the coalescence, the binary produces a steady luminosity 
Lq (which evolves only on long time-scales ^ ^orb), with roughly periodic fiuctuations of 
amplitude ALq and period tvar = (1 + z)torh about this steady mean luminosity. As argued 
in the Introduction, periodic variations could be reasonably expected if the luminosity is 
tied to the mass accretion rate, with the latter modulated on the orbital period. Even 
in the absence of such modulations, the emission could vary owing to the orbital motion 
and emission geometry of the binary (Kocsis & Loeb, in preparation). In the absence of a 
quantitative model for the electromagnetic emission, we will assume the amplitude ALq is 
unknown, and below we will ask whether a particular assumed ALq may be detectable. Since 
the residence time tres decreases continuously as the binary separations shrinks, variability 
with decreasing periods torb would be exhibited by a diminishing fraction ~ trcs/tq of bright 
quasars. 

Our main point is that an observational survey can attempt to identify such periodically 
variable sources. The total number of such periodic sources will be A^ar ~ ^res^mg, where 
Amg represents the merger rate between BHs within the survey volume (or more precisely. 



9t 



For example, the model bv lCiotti fc Ostrikeij (j200l[ ) for episodic quasar activity involves a single SMBH, 



and does not require the presence of a binary. 
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the activation rate of SMBH coalescence events). In general, the merger rate depends on 
redshift and on both BH masses, or iVjng = iV'mg(z, M, g), and should include only those 
sources with a luminosity above the survey detection threshold. To account for the latter 
condition, the light-curve of each SMBHB, Lq = Lqlt, M,q) needs to be known (here t 
could, for example, refer to the look-back time before merger). 

The merger rate N^g{z,M,q) can be modeled using the dark matter halo merger rate 
with a recipe of associating BHs with halos, and the light-curve L oft, M, q) can then be 



constrained by matching the observed quasar luminosity function (e.g. lKauffman fc Haehnelt 
20001 ). H owever, a large range of such BH population models can fit the observational 



data (e.g. iMenou. Haiman fc Narayanan! l200ll : iLippai. Frei &: Haimaru l2008bl ) . To proceed 



we instead make the simple assumption that each BH binary produces a constant mean 
luminosity of Lq = /Edd-^Edd for a total duration tq during its lifetime, where /sdd is a 
constant of ord er unity. Reasonable fiducial values appropriate to the bright quas ar phase 
are/Edd ~ 0.3 (IKoUmeier et al.ll2006l ) and, as mentioned above, ~ few x 10^ yr (iMartini 
20041 ). Note, in particular, that the quasar lifetime tq is known to be much shorter than the 
Hubble time, and N^g, which is likely determined by the galaxy merger rate, and proceeds 
on a cosmological time scale, can reasonably be assumed to be constant during tq. 

Under the above assumptions, the fraction /var of objects with luminosity Lq that 
display periodic variability on the time-scale tvar is simply given by the ratio /var = ^res/^Q- 
This ratio can be read off directly from Figures [T}171 Note that this conclusion still holds if 
the quasar emission is intermittent; we require only that the quasar is "on" for the duration 
tres when the binary orbital timescale is torb- Most importantly, under these assumptions, 
the predicted number iVvar = (tres/^Q)-?Vtot is a fixed fraction of the total number A^tot of 
quasars, and is independent of the merger rate, as long as the latter is constant during tq. 
We can then associate A^tot with the observed number of bright AGN. In particular, in the 
GW-dominated regime, we have the simple prediction 

8/3 



orb 



A^. 



tot 



10 'yr 
to 



50.2(1 + z) week 



(50) 



Note that in this equation, tvar = (1 + z)tohs is the variability time-scale as observed on 
Earth (assumed to equal the redshifted orbital time); the quasar lifetime tq is evaluated in 
the quasar's rest-frame. 

Before we proceed, we emphasize that there are many complications over the above, 
simplified picture. First, luminous quasar activity requires a near-Eddington mass-accretion 
rate, with the gas reaching within several Schwarzschild radii of one or both BHs. It is unclear 
whether abundant gas will indeed be present this close to the BHs, especially since during the 
late stages of the merger, the gas is evacuated from the inner disk by the binary's torques. 
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and the exterior gas disk is eventually unable to follow the rapidly decaying BH binary. 
Furthermore, in the final, GW-dominated regime, the tres oc scaling strictly holds only if 
any residual circumbinary gas has negligible impact on the orbital decay. This requirement 
could, in fact, contradict the assumption that the binary is producing bright emission during 
this stage. Second, in order for the emission to be periodically variable, the gas has to 
respond rapidly to the gravitational perturbations from the binary. The time-scale for this 
response is of order the local orbital time; variability on the orbital time-scale of the binary 
itself therefore again requires gas close to the binary's orbital radius. 

If the central cavity were indeed truly empty, no gas would reach the SMBHBs, and 
bright emission could not be produced. On the other hand, an empty cavity is certainly 
an idealization, and detailed models for the joint disk + binary evolution are required to 
assess the plausibility of our assumptions. Conversely, the observations envisioned here will 
constrain such models (which, again, is the main point of the present paper). 

In support of our assumptions, we note, however, that gas could be present near the BHs 
in the case of unequal masses (so that the torques are reduced), or if the disk remains thick, 
making it difficult for the binary to open and maintain a nearly e mpty central cavity. Numeri 



cal simulations indeed suggest residual gas infiow into the cavitv (lArtymowicz fc Lubowlll996 



MacFadyen fc Milosavljevidl2008l : iHayasaki et al.ll2007l : ICuadra et al.ll2009l ) , which may plau- 
sibly accrete onto the BHs ( with both BHs possibly forming their smaller individual accretion 
disks; iHayasaki et al.ll2008l ). producing non-negligible EM emission. Simulations have also 
shown, in the context of proto-planetary disks, that when the circumbinary di sk is sufficiently 
thick, the mass flow rate across the gap is increased (IDobbs-Dixon et al.ll2007l ). Such residual 
inflow onto a SMB H binary has been invoked to explain the ~12-yr periodic emission from 
the quasar OJ287 (lArtymowicz fc Lubowl Il996l ). More recent ly, the large velocity offsets 
seen in the spectrum of the quasar SDSS J092712.65+294344.0 JXomossa. Zhou fc Lul[2008h 
have be en interpreted with a similar model, inclu ding gas inflow onto a luminous SMBH 
binary (iBogdanovic. Eracleous fc Sigurdssonl 120091); a similar interpretation was invoked for 
the binary quasar candidate recently identifled by lBoroson fc Lauerl (120091 ). 



There are additional caveats that will hamper the identiflcation of any periodic sources, 
even if they exist and produce bright enough luminosity to be detectable. T he Eddington ra- 
tio of bright AGN is already known to have a signiflcant scatter (~ 0.3 dex; iKoUmeier et al. 



20061 ) . The light-curve of the merging binary is also like ly to evolve, rather than having a 



simp le "tophat" shape. It is possible, in particular (e.g. iBarger et al.ll2005l : iHopkins et al. 
20051 ) that merging SMBHs spend a signiflcantly longer time (~ 10^ yr) at lower luminosi- 
ties, (/Edd ^ !)• This will complicate the interpretation of any observed variability (i.e., 
converting the observed ratio N^av/Ntot at tvar to tres will require knowing the probability 
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distribution of Eddington ratios). This, however, can be alleviated by considering only the 
relative abundance of periodically variable objects at different values of tvar, instead of the 
absolute number of sources that show periodic variability. In this case, the only assumption 
required is that /Edd does not evolve significantly during the observed range of tvar ^ this 
should be reasonable over a factor of a ~ few range in orbital radius or in torb- Furthermore, 
even if there is a range of different BH masses, among sources with a similar luminosity, 
producing variability with the same period. Figures [1] and [2] show that more massive BHBs 
will move much more quickly through a fixed torb- Given that there are most likely fewer 
of the more massive BHBs to begin with, the set of all sources with the same tvar will be 
heavily dominated by the lowest-mass BHBs, caught at their relevant orbital radius. This 
still leaves the caveat, however, that the source is significantly sub-Eddington during the 
late stages of coalescence. In this case, the periodic sources will be harder to detect both 
because they are fainter, and also because they will also be rarer (among the long-lived and 
therefore more numerous, near-Eddington quasars with a similar luminosity). 

Another caveat is that at fixed torb and M, the distribution of q is unknown, and can 
depend on M. However, bright AGN activity is thought to be activated only in relatively 
major mergers. A smaller satellite galaxy, falling onto a larger central galaxy that is more 
than ^ 10 times more massive, may not experience the torques needed to bring its gaseous 
nucleus, with the low-mass BH, clo se to the center of the larger galaxy, for the BH-BH merger 
to take place (IHopkins et al.ll2006r): the dynamical fricti on time for small galaxies themselves 



can also be too long (e.g. iKauffman fc Haehneltl l2000l ) . and/or the smal l satellites can be 



tidally stripped before reaching the central regions of the larger galaxy (e.g. lKazantzidis et al. 



20051 ). These arguments, coupled with the well-e stablished correlations between the mass of 
a SMBH and its host galaxy (e.g. lFerraresell2002l . see also the Introduction), suggest that the 
g-distribution among binaries associated with quasars may not extend to values significantly 
below g ~ 0.1. 

Finally, for simplicity, in our estimates we have assumed circular orbits, both for 
the binary and the disk gas. It has been shown that the binary -disk interaction could 



drive both the SMBHs and the gas to have significant eccentricities (lArmitage fc Nataraian 



2005 



MacFadven fc Milosavlievidl2008l : iDotti et al.ll2008l : iHavasaki et al. 



2008 



Cuadra et al. 



20091 ). Such eccentricities should leave c haracteristic asymmet ric signatures in the modu- 
lated mass accretion rate (see Figure 8 in lHayasaki et al.l 120071 ). The resulting light-curves 
may exhibit corresponding features, which could be resolved, given sufficient time-sampling. 
In practice, allowing for eccentricities will most likely further complicate the interpretation 
of any observed period distribution, especially if the time-sampling is too coarse to explicitly 
reveal any asymmetric features. 
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3.2. Requirements of a Variability Search 

Despite the caveats hsted in the previous section, it is plausible that the periodic sources 
envisioned here exist, and we propose that they can be looked for, in a suitably designed 
survey. Most importantly. Figures [1] and [2] show that the expected variability timescale can 
be in a suitable range for a statistical detection, with a duty-cycle of t^es ^ 10^ yr over the 
range from t^ar ~day to ~yr. This suggests that such periodic sources may not be too rare. 

What will be the practical limitations for discovering the population of periodic sources? 
Clearly, there has to be a sufficient number of sources, observed over a range of variability 
time-scales for a representative statistical sampling, and the brightness variations of these 
sources must be detectable. In addition, the individual light-curves have to be sampled well 
enough to conffim their periodic nature: this will be necessary to distinguish the coalescing 
SMBH binaries from other types of variable objects. Besides discovering the periodic sources, 
the idea proposed here is to measure the dependence of iVvar on tvar - possibly to use the 
iVvar oc tvir scaliug to demonstrate that the periodic variability comes from perturbations by 
the orbital motion during the GW inspiral. For this, the survey also needs to cover at least 
a factor of several range in tvar- 

The above issues will place requirements on (i) the sensitivity and (ii) solid angle, as 
well as on the (iii) total duration and (iv) sampling rate for a survey. We can use the 
simple disk models and the idealized picture discussed above, to roughly delineate these 
requirements. For simplicity of discussion, let us assume that all sources are at z = 2. In 
reality, quasars (and therefore major BH mergers) have a broad distribution with a peak 
around this redshift; clearly this will have to be taken into account in designing an actual 
survey. For simplicity, let us also fix the mass ratio q = 1. In reality, there should be a 
distribution of v alues, perhaps in the range 0.1 < g < 1, for the mergers that activate bright 



quasar activity (IHopkins et al.ll2006l ). This would not significantly affect our conclusions. 



unless q frequently extends well below 0.1. 

Imagine a survey with a sensitivity that corresponds to detecting the periodic variability 
of BHBs with a mass Mmin at z = 2, covering a solid angle AQ. (A real survey, of course, 
will have a completeness for variability detection that is not a step function). Let us assume 
that the variable flux corresponds to a fraction r/^ar of the steady mean luminosity, ALq = 
VvsltLq = "/^var/Edd-^Edd- If the survey volume contains a total of A^tot SMBHBs with the 
luminosity Lq, then the periodic variable fraction, tres/^Q; can be determined down to the 
smallest value ~ N^l (i.e. to flnd at least one periodic source). Fixing the values of r^var/Edd 
and tq (as well as Mmin, z and q), this corresponds to a minimum variability time-scale 
tvar.min that cau be probed. Let us deflne the requirement that this minimum is tvar.min < 20 
weeks. Assuming that the longest variability time-scale of interest is around 
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year (so that the periodic nature of the variations can be convincingly demonstrated over 
a muhi-year survey), this will offer a factor of three range in tvar for mapping out the A'var 
vs. tvar dependence. For example, with the steepest possible (pure GW-driven) scaling 
iVvax oc tvar, a survey volume containing a single source with t^ar = 20 weeks would then 
contain 3*/^ ~ 20 sources with a similar luminosity but with a tvar = 60 week period. 

To fix some numbers, let us set f^dd = 0.3, ?7var = 0.1, and tq = 10^ yr. For reference, the 
Eddington luminosity of a 3 x 10^ Mq BH at z = 2, assuming a ~ 10% bolometric correction, 
corresponds to an optical magnitude of ~24 mag (in the i band). Let us also impose the 
(somewhat ad-hoc) requirement that the survey volume should contain at least A^var > 100 
sources with a detectable flux variations at the period of 60 weeks. In the GW-driven stage, 
there will then be at least 5 detectable periodic sources with a period of < 20 weeks; in the 
gas-driven regime, where the scaling /^ar vs. tvar is flatter, there will be a larger number of 
20-week period sources. 

In Figure [H the curves show the sky coverage required to satisfy these criteria, as 
a function of the i-band variable magnitude corresponding to the detection limit of the 
survey. The BH masses producing the corresponding steady i magnitude (which, in our 
fiducial model, is 2.5mag brighter than the variable magnitu de) are shown o n the t op axis. 
This figure assumes q = l- We used the fitting formula by iHopkins et al.l (l2007bl ) for the 
bolometric quasar luminosity function (LF) dcj) / dL{z , L) to compute the the total number 
A^tot of quasars at z = 2, per solid angle Afi, in a redshift range of Az = 1, i.e. A^tot = 
{^z/S.VL){dPV / dzdVL) {d(f)/dL)dL, where (d'^V/dzdfl) is the cosmological volume element, 
and Lmin is the bolometric luminosity corresponding to the steady magnitude threshold i. 
We then used equation (l50l) for /var to obtain the total number iVvar = /varA^tot of variable 
sources at observed period of tvar = 20 weeks and at t^ar = 60 weeks. Requiring A^var(^var = 60 
weeks) > 100 then yields the solid angle AQ as a function of i. Note that the quasar LF 
is almost a pure power- law up to Mbh ^ 10^ Mq. The break between z = 26 - 27 mag m 
the solid curves corresponds to the transition between GW and gas-driven orbital decay. 
In particular, the figure shows that SMBHs with a mass above/below ~ lO^Mo are in the 
G W/gas-driven regime, respectively. 

Figure [H] shows that there is a clear trade-off between survey depth and area: the 
required sky coverage scales with the survey flux limit approximately as AQ oc with a 
steepening for shallow surveys with limiting magnitudes i < 22.5 (due to the decline at the 
bright end of the quasar LF), and a flattening for very deep surveys with limiting magnitudes 
i > 26.5 (because the SMBHBs are in the gas-driven regime and their residence times at 
fixed torb are shorter than in the pure GW-driven regime). 

From Figure [HI we conclude that, for example, a 1 sq. degree survey, detecting SMBHBs 
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whose steady luminosity is i = 23.3 mag, with a variabihty at the level oii = 25.8 mag, with 
sufficient sampling and duration to cover periods of 20-60 weeks, represents an example for 
the minimum specification for the survey parameters (i)-(iv). In this example, the mass of 
the BHs being detected is ~ 2 x 10'' Mq. Figures [T] and [2] show that at orbital periods of 
(20 — 60)/(l + 2;) = 7 — 20 weeks, these BHs are all in the GW-driven regime when g ~ 1, 
but may be in the viscosity-driven regime for g -C 1. Surveys that go deeper and cover a 
smaller area will begin probing the gas-driven evolutionary stages. 

In Figure [9l we examine how the required sky coverage changes when the parameters q, 
Vvar or /Edd are modified. The middle (black) curve shows the sky coverage required to find 
20 sources at tvar = 35 weeks (intermediate between the red and black curves in Figure [8]), 
with our fiducial parameters, q = I, //var = 0.1, /sdd = 0.3. The green curve corresponds 
to changing the mass ratio to g = 0.01; this increases/decreases the residence time in the 
GW/gas-driven regimes relative to the q = 1 case (compare Figs. [T]and[2]), and therefore 
reduces/increases the required solid angle coverage. The top pair of (blue) curves in Figure [9] 
show variations when either r^var or /sdd is decreased by a factor of 10 (to 0.03 or 0.01, upper 
and lower of the pair, respectively). Note that the survey volume requirement is more 
sensitive to r^^ar (whereas the critical BH mass is equally sensitive to either). Similarly, the 
bottom pair of (red) curves show variations when either r/var or f^dd is increased by a factor 
of 10. The small break visible at ~ 25.5 mag in this case corresponds to the transition from 
the middle to the outer disk region for ^ 2 x 10^ M© SMBHBs at torb = 35/(1 + 2) ^ 11.5 
weeks (see Fig. [T]). Each dashed curve shows the mass of the SMBHB corresponding to the 
effective i magnitude limit (BH masses are labeled on the right y-axis). The required survey 
volume also shifts linearly with the assumed total quasar lifetime tq. 

If we knew where to look (i.e., if LISA delivers a candidate for an on-going merger, 
with sufficiently accurate localization on the sky), it would be possible to perform a deep, 
targeted observation for variability on short time-scales; between several minut es up to ~ 10 



hours within the last ~ month of merger (this possibility is discussed in detail in lKocsis et al. 



20081 ). However, each source will spend only a ~month at such short variability time-scales. 



and a random search, in the absence of a preferred direction on the sky, would then have to 
monitor > 10'^(tQ/10^yr) AGN to find a single example of such a late-stage periodic source. 
Alternatively, one may monitor ~ 10^(tQ/10^yr) AGN for ~ 10 years, to look for (slowly 
evolving) periods, on the timescale of ~ a day. The slow decrease in the period, which would 
be a smoking gun for GW-inspiral, will be challenging to observe in real time for individual 
objects. 
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3.3. Constraints from Existing and Future Data 



Existing observations from radio to X-ray bands liave shown that the luminosity of 
quasars and other active galactic nuclei vari es on time-scales from h ours to sev eral years 
(see, e.g., the articles in iGaskell et al.l l2006l or the recent review by iFanI |2005| ). In fact, 
variability often aids in the identification of AGN (and may conversely be a major obstacle 
in identifying the periodic signal proposed here). While variability is detected in a large 
fraction of all AGN, there are only a handful of sources whose str ucture functi on shows 
clear periodic variability on long (> weeks) time-scales (see, e.g. iRiegerl 120071 for a re- 
view focusing on searches for periodic variability). Examples include a handful of blazars, 
whose historical light-cur ves show periodic ou t bursts on timescales o f a year to a (decade 



or even longer (see . e.g.. ISillanpaa et al 



Raiteri et al.l l200ll : iQian et al. 



2007 



1988 



Tao et al. 



2008 



and references in these papers, for 
individual objects). iNeugebauer fc Matthews! (119991 ) monitored 25 low-redshift, optically 



iu. Xie fc Bail Il995l : iLiu. Liu fc Xid Il997 



selected quasars for variability over several decades in infrared bands. They identified one 
object, the radio-loud quasar PG 1535+547, whose structure function shows a periodic com- 
ponent with a period of ~ 10 yr. This source has a bolometric luminosity of Lboi = 10^^ '^'^Lq, 
implying a BH mass of ~ 3 x 10^ (/Edd/O.3)~^M0. Figure [1] shows that equal-mass SMBH 
binaries with this mass are in the GW-driven regime, and /var ~ 20% may exhibit a 10-year 
period. In comparison. Figure [2] shows that unequal-mass binaries with this total mass may 
be in the gas-driven regime, and periodic variability may be exceedingly rare. Thus, we 
conclude that the identification of one periodic object is roughly consistent with it being an 
example of a GW- driven, near-equal mass b i nary. However, there are only 4 objects in the 
sample studied by INeugebauer &: Matthews! (Il999l ) with luminosities above Lboi = lO^^Lp 
prohibiting robust conclusions. 



0) 



Very large area variability surveys, such as in the S loan Digital Sky Survey ( SPSS), are 
shallow, and detect variability down to only i ^ 20 mag (IVanden Berk et al.!l2004l ). With our 
fiducial /var = 0.1, this corresponds to periodic AGN whose steady luminosities are i ~ 17.5 
mag. Figure [S] shows that such a survey would have missed the periodic variations discussed 
above. Figure [9] shows that if the variable fraction is much larger, r/var ^ 0.3, then ~ 100 
variable sources with a yearly period would be detectable - but their periodic nature could 
be demonstrated only with sufficient time-sampling, extending over a decade. 

A recent, much deeper optical survey by the Subaru telescope ( Morokuma et al.ll2008a| jbl) 
for variable objects provides interesting constraints on the scenario envisioned here. The 
completeness function in this survey, defined as the probability to detect fiux variations of 
an object wit h a variable component i , goes from ~unity to ~zero between i ^ 25 to i ^ 26 
(see Fig. 8 in lMorokuma et al.!l2008a! ). In the fiducial case with f^dd = 0.3 and r^^ar = 0.1, 
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25.5 mag corresponds to the mean steady magnitude 
= 2.5 X 10^ Mq. At z = 2, the Subaru s urvey has a 



the hmiting variabihty magnitude i = 
of i = 23 mag, and BH mass of M 
completeness of 0.5 at this magnitude (see Fig. 11 in iMorokuma et al.l l2008al ) and covers 
an area of 0.9 sq. degrees. This combin ation of sensitivity ari d area hes very close (just 
below) the curves in Figure [81 Using the iHopkins et al.l (l2007bl ) quasar LF, and assuming 
a completeness of 0.5, we find that the Subaru survey should detect 440 AGN; this is in 
nearly exact agreement with their quoted result (489 deg"^). We further find that of these 
sources, ~ 0.6, 12, 61 would vary with observed periods of 20 weeks, 60 weeks, and 1000 



days ( adopting a pro bability of 0.5 for detecting variability, from Fig. 8 in iMorokuma et al. 



2008al ). Figure 12 in IMorokuma et al.l (l2008bl ) shows that they found several dozen sources 
that varied, at least once in their life, on all of these timescales. Unfortunately, we do not 
know whether these sources are periodic or not, and therefore the Subaru survey results 
represent only an upper limit on the fraction of periodic sources. Nevertheless, this already 
suggests that the 2.5 x IO^Mq BHs at the limit of the survey can not produce variability at 
the level significantly exceeding our fiducial 0.03LEdd- 

The ultra-deep Hubbl e Space Telescope (HST) variability surveys (see, e.g., the recent 
review by ISarajedinil |2008| . and references therein) discovered galaxies whose nuclei varied 
by magnitudes down ioV^, 27.5. The observations were taken a year apart in the Hubble 
Deep Field North (HDFN) and the Groth Survey Strip (GSS), whose areas are ~ 3 x 10~^ 
and ~ 0.1 sq. degrees, respectively. While the solid angle of the HDFN dataset is too small 
to yield useful constraints, the fiducial case with /Edd = 0.3 and r/var = 0.1 in Figure [S] shows 
that the GSS dataset just reaches the sensitivity/area combination of ~ 27 mag and ~ 0.1 sq. 
degrees required to find flux variations from Mbh ~ lO^M© SMBHBs. Appro ximately 4.5% of 
AGN were found to vary by magnitudes down to V ~ 27 in this dataset (see lSarajedini et al. 



20061 . for more details), suggesting that < 5% of AGN containing SMBHBs with this mass 
can produce variability at the ~ 0.03LEdd level. 

AGN are also known to vary on long times-scales in X-ray bands. Systematic and unbi- 
ased variability surveys sensitive to times-scales of we eks, such as those in soft X-rays in the 
ROSAT all sky survey (iFuhrmeister fc Schmittll2003l ) or in hard X-rays in Swift/BAT data 
( iBeckmann et al.ll2007l ) howev er, have been restricted to th e brightest AGN, while deeper sur- 
veys, such as those by RX TE (IMarkowitz fc Edelson 1200 ll) , of the Chandra Deep Field North 



(iBauer et al.ll2003l . |200J) and South (IPaolillo et al. 



2004h . and by XMM flPapadakis et al. 



20081 ) have only monitored up to a few hundred sources. These observations do suggest that 
a large fraction of AGNs vary in X-ray bands on time scales of a day to a year, but whether 
the variations are periodic h ave not been determined . In the 9-month duration observations 
covered by Swift/BAT data, IBeckmann et al.l (120071 ) find a strong anti-correlation between 



luminosity and variability (with no source with luminosity Lx > 5 x 10 erg s showing 
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significant variability); iPapadakis et al.l (120081 ) report a similar trend from an X-ray vari- 
ability analysis of 66 AGN in the Lockman Hole. These findings would be consistent with 
the trend that the most massive SMBHBs (> 10^ Mq) spend less time at a fixed orbita l 
timescale of torb ~ 10 weeks (see Figures [1] and [2]) . The results of iBeckmann et al.l (120071 ) 
suggest that absorbed sources vary more than unabsorbed ones, which may be particularly 
relevant for finding the periodic SMBH binary sources envisioned here, which are undergoing 
the last stages of their merger, and may be heavily obscured and visible primarily in X-ray 
bands. 

While the deep existing optical surveys come close to placing useful constraints on the 
scenario envisioned here, future surveys, designed to uncover source populations with pe- 
riodic variations on times-scales of tens of weeks, should be able to either discover these 
populations, or place stringent limits on their existence. Many large optical/IR surveys are 
being planned or buil t, motivated l argely by finding type la supernovae (SNe) for cosmolog- 



ical studies (see, e.g.. IStubbs 



2008. for a recent review). The most ambitious of these, such 



as LSST and Pan-STARRS-4lj will be all-sky surveys, and should be able detect variability 
to 26 — 27mag, allowing detectio ns well beyond the mo st pessimistic case shown in Figure [9l 
The proposed ALPACA survey (ICorasaniti et al.ll2006l ) will cover 1,000 sq. degrees to 23-25 
mag in 5 optical bands, and would already reach the sensitivity /area combination probing 
these pessimistic scenarios. 



3.4. Other Detection Methods 



In addition to producing periodic variability, there could be several other methods to 
prove or disprove the presence of a SMBH binary. First, the orbital motion of the binary may 
cause relative shifts in the quasar's emission lines. For example, in a configuration in which 
the broad lines arise from gas close to one of the two (moving) BHs, and the narrow lines arise 
from material farther away, which is close to rest at th e systemic redshift, such a shift could 
arise between the narrow and broad emission lines (iBogdanovic. Eracleous fc Sigurdsson 
20091 ) . Similarly, if both BHs ca rry their own accretion disks, extending to a few Schwarzschild 
radii, and produce broad lines (IHayasaki et al.ll2007l ). then there could be two sets of broad 
emission lines, super-imposed with a similar relative velocity shift. The magnitude of these 
shifts may be of order the orbital velocity {v ~ 6, 000 km/s at ~ 1, 000 Schwarzschild radii), 
which could be detectable either in individual objects, or else statistically for the population. 



see www.lsst.org and www.pan-starrs.ifa.hawaii.edu 
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Boroson fc Lauerl (120091 ) recently reported a candidate SMBH binary, with two sets 
of broad emission lines separated by 3, 500km s~^. The spectrum of this s ource can als o 
be interpreted with a single BH+disk system (IHalpern &: Filippenkd Il988l : iGaskelll l2009l ): 
indeed, this interpretati on is favored by the lack of any change in the velocity offset over 
the course of ~ 1 year (jChornock et al.l 120091 ) . Nevertheless, it is interesting to note that, 
with the binary parameters reported for this source (assuming random orientation, and an 
expected orbital speed of ^ 6, 000km s'^), M ^ Mi = lO^-^ Mq, Ma = lO^-^ Mq, q = 0.025, 
torh ~ 100 years, and R ~ lO^Rs, we find the evolutionary track of the proposed system to 
be virtually indistinguishable from the M = 10^ M©, q = 0.01 case shown in Figure |H At its 
currently observed orbital separation of ~ lO^Rs, the binary would be in the gas-driven 
regime, close to outer radius of the formally gravitationally stable disk (i.e., the system is 
just outside the marked Q = 1 point in this figure), with a residence time of ~ 10^ years. 
This could indeed make this observed separation common among quasars with M = 10^ Mq 
SMBHBs. The figure also shows that the residence time at fixed orbital velocity decreases 
steeply with BH mass, suggesting that fainter binary quasars with similar orbital speeds 
would be much less common. 



In the last stages of coalescence, the GWs emitted by such a system would induce 
periodic modulations in the arrival times of pulses from background radio pulsars; at 200ns 
timing sensiti vity, these modul ations would be detectable from SMBHBs out to a distance 
of ~ 20Mpc (jjenet et al.ll2004l . note that this study already applies the idea to the source 
3C66B mentioned in § 1, whose elliptical motion was interpreted as due to a SMBHBs, and 
rules out the SMBHB hypothesis). 



4. Summary and Conclusions 

In this paper, we followed the evolution of SMBH binaries, starting from large separa- 
tions, to coalescence. We find, in agreement with earlier works, that the orbital decay is 
initially generically driven by viscous binary-disk interactions, whereas GWs dominate the 
last stages. In a refinement of earlier results, we also find that just prior to the transition to 
GW-driven evolution, the viscous orbital decay is generically in the "secondary-dominated" 
Type II migration regime (the mass of the secondary is larger than the enclosed disk mass). 
This is slower than the disk-dominated Type II migration that has sometimes been assumed 
in the past, and, as a result, SMBH binaries spend a significant fraction of their time at 
orbital periods of ~days to ~ a year, where they may not be rare, and may be identifiable. 
We emphasized the large uncertainties in the residence times in this regime - for example, 
time-dependent disk models predict even slower decay. We also find that observations of 
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BHs with a mass range of 10^ — 10^ Mq over this range of periods could find binaries lo- 
cated in all three physically distinct regions of the circumbinary disk. Thus, several aspect 
of disk physics could potentially also be probed in future observations of a population of 
SMBH binaries. We also find that viscous processes may contribute to the orbital decay rate 
even after the binaries enter LISA^s frequency range, for low- and/or unequal-mass binaries 
(M < 10^ M0 or g ^ 0.01). While viscous processes are strongly sub-dominant for rapidly 
evolving "inspiral" sources, detected during the last few years of their coalescence, the pres- 
ence of the gaseous disk could reduce any background of unresolved stationary sources at 
frequencies near the low-frequency end of the LISA range (~ 10~^ mHz) . 

We considered the possibility that there may be a one-to-one correspondence between 
the activation of luminous AGN and SMBH coalescences, with a fraction of AGN exhibiting 
periodic flux variations. Given that the interpretation of individual SMBHB candidates have 
so far remained ambiguous, we proposed that a statistically large sample should aid in the 
identification of these binary BH sources. Our main conclusion is that future surveys in 
optical and X-ray bands, which can be sensitive to periodic variations in the emission from 
~ 10^ — 10^ Mq supermassive black hole binaries, on timescales of t^ar ~ tens of weeks, at 
the level of ~ 1% of the Eddington luminosity, could look for a population of such sources, 
with the aim of determining the fraction /^ar of sources, at a given redshift and luminosity, 
as a function of tyar- In our simplified models for the binary-disk interaction, this time-scale 
of tens of weeks corresponds to the orbital time when binaries with M ~ 10^ Mq make 
their transition from viscous to GW-driven evolution. In the latter regime, for sources with 
M ^ 10^ Mq, gravitational radiation predicts the scahng /^ar oc tl^. The discovery of a 
population of periodic sources whose abundance obeys this scaling would confirm that the 
orbital decay is indeed driven by GWs, and also that circumbinary gas is present at small 
orbital radii and is being perturbed by the BHs. Deviations from the tvar power-law for 
lower-mass BHs would constrain the structure of the circumbinary gas disk and viscosity- 
driven orbital decay. 

There is certainly a possibility that the periodic sources envisioned here do not exist (e.g., 
because the SMBH binary does not produce bright and variable emission during its GW- 
emitting stage, at orbital separations of ~ lO'^ Schwarzschild radii). Nevertheless, we argued 
that existing surveys already approach the required combination of sky coverage and depth, 
and future surveys, designed to make observations for several years, with a sampling rate of 
a few days, could yield a positive detection and identify periodic source populations. This 
would bring rich scientific rewards, possibly including the indirect detection of gravitational 
waves, driving the orbital decay of these sources. 
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Fig. 1. — The evolution of equal-mass (g = 1) SMBH binaries, embedded in a steady 
circumbinary disk, from large to small orbital separations. The figure shows the residence 
time tres = —R/ (dR/dt) that each binary spends at the radius where the orbital time is torb- 
The four curves correspond to binaries with total masses of M = 10^, 10^, 10^, and 10^ M©, 
as labeled. The large dots denote the critical radius beyond which the assumed circumbinary 
disk is unstable to fragmentation (Toomre parameter Q < 1). Similarly, triangles denote radii 
beyond which the disk may be susceptible to ionization instabilities (the gas temperature falls 
below lO^K). In each case, blue/red colors indicate whether the disk mass enclosed within the 
binary's orbit is larger /smaller than that mass of the secondary. The dotted/dashed/solid 
portion of each curve indicates the outer /middle/inner disk region, respectively (as defined 
in § 12.21) . For a binary located at redshift z, the redshifted values of tj-es and torb (as measured 
on Earth), should be multiphed by a factor of {1 + z). 
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Fig. 2. — The figure shows the residence time t^es as in Figure [H but for unequal-mass 
binaries (g = 0.01). 
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Fig. 3. — The residence time for equal-mass binaries, as in Figure [3l except that tres is here 
shown as a function of orbital separation R, in units of the Schwarzschild radius Rs- For 
reference, the x axis labels on the top show the orbital velocity corresponding to each value 
of R/Rs. 
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Fig. 4. — The residence time tres as in Figure[3l except for unequal-mass binaries (g = 0.01). 
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Fig. 5. — The residence time for equal-mass binaries, as in Figure HI except that for each 
binary, we also show the results for 6=1 (top curves), in addition to the 6 = case (bottom 
curves, reproduced from Fig. [3]). In the 6=1 case, the viscosity is proportional to the gas 
pressure, rather than the total pressure. Once the binary approaches the radiation-pressure 
dominated regime (shown in solid), the evolution is therefore slower than in the 6 = case. 
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Fig. 6. — The residence time for equal-mass binaries, as in Figure [Tj but with additional 
curves showing the evolution expected in a scenario with a time-dependent disk. The dot- 
ted rn agenta curves were calculated based on the model by llvanov. Papaloizou. fc Polnarev 
(]l999l ). assuming that the binary-disk interaction turns on when the disk dominance pa- 
rameter reaches qb = I (prior to this, the steady-disk solution is applied). As expected, the 
pile-up of material in the time-dependent disk slows down the decay of the binary. 
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Fig. 7. — The residence time tr^s as in Figure [6 , exce pt that the binary-disk interaction, 
as modeled by llvanov. Papaloizou. fc PolnarevI (Il999l ). is turned on later, when the disk 
dominance parameter drops below qs = 0.01. Prior to this, the disk-dominated, steady- 
disk solution is applied. 
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Fig. 8. — The sky coverage required to find a population of periodic sources, assuming our 
fiducial set of source parameters (see text). The upper (black) curve shows the sohd angle 
required to find at least ~ 100 periodic sources between 1.5 < z < 2.5 with an observed period 
of 60 weeks, as a function of the limit of the survey for the variable i-band luminosity. The 
lower (red) curve shows the solid angle required for at least 5 periodic sources with t^ar = 20 
weeks. In the GW-driven regime, these two criteria coincide. The upper labels show the 
mass of the SMBHB producing the corresponding steady i magnitude (assumed here to be 
2.5 mag brighter than the variable magnitude). As shown by the break between 26 - 27 mag 
in the curves, BHs with a mass above/below ~ 10 ''Mq are in the GW/gas-driven regime, 
respectively. 



-45- 



lOg[(^varfEdd/OO3)M5h/M0] 



9 8 7 6 




i mag 



Fig. 9. — This figure shows how the sky coverage required to find a population of periodic 
sources depends on variations in our fiducial set of source parameters. The middle (black) 
curve shows the sky coverage required to find 20 sources at tvar = 35 weeks, with the same 
set of fiducial parameters, q = I, 77var = 0.1, /Edd = 0.3, as in Figure [Hi The green curve 
corresponds to changing the mass ratio to g = 0.01; the top/bottom (blue/red) pair of 
curves show variations when either r^var or f^dd is decreased/increased by a factor of 10. 
Note that the survey volume requirement is more sensitive to r^^ar (whereas the mass of the 
smallest detectable periodically variable SMBHB, shown on the top axis, is equally sensitive 
to either). 
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